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Abstract Equation of motion for the galactic tide is derived under the assumption 
of cylindrically symmetric gravitational potential of the Galaxy. The paper considers 
galactic tide both for the galactic plane x— and y— components and also for the normal 
z— component. The x~ and y— components of the acceleration come not only from 
the X— and y— components of the position of a body, but also from its 2:— component 
of the position vector. 

Values of the Oort constants are A = (14.2 ± 0.5) km kpc^^ and B — {— 
12.4 ± 0.5) km s^^ kpc~^. Mass density in the solar neighborhood, 30 pc above 
the galactic equatorial plane, equals to (0.117 ± 0.005) Mq pc~^ . The result for the 
acceleration is written in the form easily applicable to Solar System studies, to the 
evolution of comets in the Oort cloud. 
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1 Introduction 

Global galactic gravitational field influences relative motion of two close bodies in 
the form of galactic tide. This relative motion is important in Solar System studies on 
evolution of comets within the Oort cloud of comets. This paper derives the equation of 
motion for galactic tide under the assumption of cylindrically symmetric gravitational 
potential of the Galaxy. The effect of the normal z— component is well understood 
in the case when motion of one of the bodies is fixed to galactic equatorial plane 
(see, e.g., Mihalas and McRae Routly 1968, pp. 221-222; analytic approach to the 
secular evolution of orbital elements in the gravitational field of the Sun is presented 
in Klacka and Gajdosfk 1999). Acceleration in the galactic plane comes not only from 
the X— and y— components of the galactic field, but also from the 2:— component of 
the galactic components. Similarly, acceleration in the z— component comes also from 
the X— and y— components of the galactic field acceleration/tide. The acceleration 
in the plane parallel to the galactic equatorial plane and the acceleration normal to 
the galactic equatorial plane are presented in the relevant general form for two sets of 
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galactic components; one case corresponds to Dauphole et al. (1996) potential model 
for galactic bulge, disk and halo, the other case considers potential galactic bulge of 
Dauphole et al. (1996), galactic disk is represented by mass density function (Maoz 
2007) and galactic halo is given by flat rotation curve. 



2 Values of the Oort constants 

In order to have a realistic model for galactic gravitational force in disposal, one needs 
to know also the values of Oort constants. They refer for the region of the Sun, in 
galactocentric distance about 8 kpc. 

There is relatively great uncertainty in the values of the Oort constants. The values 
with errors are given in Table 1, where also sources of the values are given. 



A 


B 


source 


[km s^^ kpc^^] 


[kro kpc^^] 




15.0 ± 0.8 


- iO.O ± 1.2 


Kulikovbkij (19«r), p. <)()) 


14.4 ± 1.2 


- 12.0 ± 2.8 


Kerr and Lynden-Bell (1986) 


11.3 ± 1.1 


- 13.9 ± 0.9 


Hanson (1987) 


14.8 ± 0.8 


- 12.4 ± 0.6 


Feast and Whitelock (1997) 


14.5 ± 1.5 


- 12 ± 3 


Barbierl (2007, p. 167) 



Table 1 Values of the Oort constants taJsen from various sources. 



The values presented by Hanson (1987) are given also by Oiling and Merrifield 
(1998), the values presented by Feast and Whitelock (1997) can be found also in Sparke 
and Gallagher (2007, pp. 92-93) and Carroll and Ostlie (2007, p. 913). 

The scatter of the values of A and B is large. The best estimates of the quantities 
A and B are the weighted averages 

/I 

Awav - jv ' 



Bwav — 



(1) 



where the weight wa i {wg 
corresponding uncertainty 
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Using the above presented values, we can calculate, on the basis of Eqs. (l)-(3), the 
most probable values of A and B (we omit the index "wav" ) : 

A = (14.2 ± 0.5) km s"^ fcpc"^ , 

B = {- 12.4 ± 0.5) km s"^ kpc~'^ . (4) 

The values are consistent with the statements that "A has a value of about 14 km kpc 
and "B takes a value around — 12 km kpc~^" (Phillipps 2005, p. 101). However, 
these values significantly difi'er from the lAU recommended values A = 15 km kpc~^. 
B = — 10 km. s^^ kpc^^ (Karttunen et al. 2007, p. 358). Also the results obtained 
from the model of Dauphole ct al. (1996) are not fully consistent with the results pre- 
sented in Eq. (4). The model of Dauphole et al. (1996) yields A = 14.2 km kpc~^ , 
B = — 13.9 km kpc^^ . Similarly, the cepheids data of Karimova and Pavlovskaja 
(1974) (see also Kulikovskij 1985, p. 98) yield rotation curve v{R) = uq = (243 ± 1) 
km s^^ , which leads to ^ = — i? = 15.2 km kpc~^ for R = 8 kpc [the angular 
velocity is ij{R) = vq/R, A{R) = {v{R)/R - dv{R)/dR) / 2 = uj{R) / 2, in this special 
case when v{R) = vq: if ijj{R) is a convex function, then A{R) is the convex function 
a;(i^) / 2]. 



3 Models of Galaxy 

Galaxy is standardly considered to consist of three components, galactic bulge, galactic 
disk and galactic halo (sec, e.g.. Sec. 6 in Maoz 2007). If one does not take into account 
spiral structure of the Galajcy, then simple models can be created. 

The model of Dauphole et al. (1996) considers spherical symmetry for galactic 
bulge and halo and cylindrical symmetry for galactic disk. Disadvantage of the model 
is that it does not yield values of Oort constants corresponding to Eqs. (4). Moreover, 
the model cannot be used as a realistic model of Galaxy for galactocentric distances 
larger than about 40 kpc, since the model produces a decreasing rotation curve for 
these distances. 

We can create a new simple model, better consistent with Eqs. (4) and with flat 
rotation curve of the Galaxy. The new model considers galactic bulge of Dauphole et 
al. (1996) and its gravitational potential is 



^bulgeir) 



Mb = 1.3955 X 10 
bf, = 0.35 kpc , 



> 



(5) 



where G is the gravitational constant. 

Galactic disk of our new model is represented by mass density function. We can 
take relevant information from, e.g., Maoz (2007, Sec. 6.1.1.1): "The Galactic disk has 
a mass distribution that falls exponentially with both distance R from the center and 
height z above or below the plane of the disk: 



exp 
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The scale length of the disk, R^t = (3.5 ± 0.5) kpc, and hence at = 8 kpc the Sun is 
in the outer regions of the Galaxy. The characteristic scale height is = 330 pc for the 
lower-mass (older) stars in the disk and = 160 pc for the gas-and-dust disk. The Sun 
is located at z = 30 pc above the midplanc of the disk. The mass of the disk within one 
scale radius is MfUsk ~ most of it is in stars (and about 10 % in gas)." Eq. 

(3) is consistent with the statement presented by Bertin (2000, p. 33): "The luminosity 
profiles of galaxy disks are approximately exponential, being reasonably well fitted by 
the law I{R) = Iq exp(— _R//i), which has two scale parameters, the central brightness 
Iq and the exponential length h." . 

Finally, the simple model of galactic halo is given by a flat rotation curve. The 
rotation curve is given by circular speed Vhaioi''') for spherical halo: 



[vhalo{r)f = ■f'ff S 1 - a ^ arctg ( ^ ) " (1 - «) exp 



r \^ 



vji = 220 km s ^ , 
a = 0.174 , 

an = 0.04383 kpc , 

bH = 37.3760 kpc , (7) 

(compare with Sparke and Gallagher 2007, p. 95). The numerical values of a, a^j and 
bn are calculated on the basis of Eqs. (4)-(7) for the given value of Vfj (although a 
different value can be used), [vhalo{R)? = ["GalaxyiR)? - [vbulge{R)f - bdisfc(-R)]^, 
VGalaxy{R) = (^ " B)R, dvcaiaxyiR) / ^R = ~ [A + B); moreover, the values yield 
minimum potential energy of the halo and also minimum radius of the Galaxy for a 
given mass of the Galaxy. 



4 Motion in GalcLxy — general description 



Let us consider an approximation when global galactic gravitational field is described by 
cylindrically symmetric potential ^{R,z), R being distance from the axis of rotation 
and z the coordinate of a body above/below the galactic equatorial plane (z = 
corresponds to the galactic equatorial plane; right-handed system x — y — z has its 
origin at the cente r of the Galaxy, z is positively oriented toward the north pole of 
the Galaxy; R = -\- j/^). The galactic gravitational potential is generated by mass 
distribution within the Galaxy. We can write 

^{R,z)=^^ulae{R.z) + ^haloiR.z) + 0disk{R,z) ■ (8) 

Acceleration of the body is given by the following equations in cartesian coordinates: 

d${R,z) X 

dR R ' 
d<l>iR, z) Y_ 
dR R ' 



(fx 


d${R, z) 


df^ 




dx 




d<l>{R, z) 


df^ 




dY 






'd${R,z)' 


dt^ 




dz 


R = 


+ y2 _ 



z=Z 



(9) 
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For the case represented by Eqs. (5)-(7) we have 



r = A/i?2 + «2 ^ 
^disk{R^ 2:) is found from the solution of Poisson's equation 

and, finally, 

r = ^Jb? + 2^ , 
d^haloir) _ [vhaloir)? 



(10) 



(11) 



dr 



(12) 



5 Motion near the galactic equator 

On the basis of Eqs. (5)-(7) and (10)-(12) we can write 



= ^ (1 + Oi 

Rq Rq 
i = bulge, disk, halo , 



z + 



and. 



Ho 

,2 



2 



, 
i 

VG = vo = v{Ro) = {A - B) Ro . 
Using numerical values for the model, we can write 



OR 



iR,z) 



{A - BY Ro 1 



-Ro 

A = 0.124 fcpc" 
= 1.586 kpc~ 
Ro = 8.0 kpc . 



Bl Z + - B2 Z 



(13) 



(14) 



(15) 



We have added also the value of i?Oi the galactocentric distance of the Sun. The values 

of A and B arc given in Eqs. (4). 

The rotation curve is given by the circular velocity v{R) = ^ R d${R, 0) / dR. Eqs. 
(9) can be written, then: d^X/di^ = -{v'^/R^)X and d^Y/dt^ = -{v'^/R^)Y for x- 
and y-componcnts in the plane z = Q. Using two close points, (Xo, Yq, Zq) and {X, Y, 
Z), X = Xo+t Y = Yo+r], Z = Zo+C, one can write {Xoi+Yon+ZoQ/ ^/R^+Z^ « 
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(Xq^ + lo^ + ^oC)/-Ro for the difTerencc between magnitudes of their position vectors, 
if higher orders in Zq/Rq, rj and C are neglected. For the galactic plane = 
VQ{l + 2{vQ/vo)(XQ^ + Yor})/Ro}, where the prime denotes differentiation with respect 
to R {vq = v{Rq), v'q = [dv{R)/dR]jig) and, again, higher orders in ^ and rj are 
neglected. We are dealing only with \z\ <C 1 kpc; Zq = 0.03 kpc, at present (Maoz 
2007), -Ro = 8 kpc. 

The total action of all galactic components can be summarized (see Eqs. 9): 



d^X 
dt^ 



d^Y 
dt^ 



4 

Rl 



Xo 



4 

Rl 



- Yo 



Xo+? + 2 I Ra^ 

n {zl + 2 zo c) 

yo + 77 + 2(7io^- 
A {zl + 2ZoC) 



r2 {zt + 4 zl c)] } 



Rl 



(4 



<)]} 



(16) 



if higher orders in ^, 77, C, are neglected. 



6 Relative acceleration 



The relative acceleration of the body with respect to the Sun, if (Xq, Yq, Zq) represents 
position vector of the Sun, is: d^^/dt'^ = - GMq^/v^ + d^ X/ dt^ - d^ Xq/ dt^ , d^rj/dt'^ 
= - GMQTi/r^ + d^Y/di^ - d^Yo/dt^, or 



dt^ 



df^ 



GMq ^ 



-2XoZq (a - r2 zl ) c} , 



Rl 



XqYq 
Rl 



GMq ^ 



Rl 



Rl 



-2YqZq (a - r^zl) , 

r = + + C2 . 

Using Oort constants A and B, fulfilling the relations A 
= -v'o, Eq. (17) yields 



Vo 



(17) 

wo = vq/Rq, A + B 



GMq 



- —^£, + {A-B)[A^-B + 2Acos{-2ujQt)]^ 
+ 2A{A - B) sin (-2 ujot) rj 

+ 2 {A- Bf (A - A ^0 ) Ro Zq cos(- woi) C , 



dt^ 



7] + 2A{A - B) sin (-2 woi) C 
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+ [A - B)[A + B -2A cos (-2 uqI)] r/ 

+ 2 {A- Bf (A - r2Zl) Ro Zo sin(- wq*) C , 

r = + + , (18) 

where the sign minus at angular velocity (— ljq) denotes negative orientation of the 
galactic rotation (clockwise orientation/direction of the solar motion with respect to 
the center of the Galaxy), and, also the relations Xq = Rq cos(— ojot) and Yq = Rq 
sin(— Loot) are used. 



6.1 0— component of the acceleration 

Eqs. (18) represent x— and y— components of acceleration. In order to be the system 
of differential equations complete, we need also the 0— component of the acceleration. 
It can be obtained from the Poisson's equation 



RdR 



+ 



= 4 71- G ^(-R,^;) 



(19) 



Taking into account that Rd<P/dR = [v{R)Y and [v{R) / R]ji^^ = luq ^ A - B, [dv{R) /di?]_R(, 
= — {A + B), where A and B are the Oort constants, and, the component of the 
acceleration is az = — d$/dz, Eq. (19) yields 



dz 



Ait G g{Ro,z) - 2 (A 
'dQ{R,z) 



B^) 



- 47r G 



dR 



Ro 



(20) 



The first term in Eq. (14) is dominant for the region around the Sun. Taking q{Ro, z) 
= q{Ro,zo) in the right-hand side of Eq. (20), one finally obtains 



[4 TT G g{Ro,ZQ) 
dg{R, zo) 



- 4nG 



dR 



- 2{A^ - B^)] z 

Kq 



-Ro 



iv-Yo)} 



(21) 



We can take g{Ro, zo) as the sum of densities coming from the disk, halo and bulge. 
The density coming firom the bulge is of negligible importance (see, e.g., the bulge 

model of Dauphole et al. 1996). The contribution of the disk from the Maoz model 
yields Qdisk = (0.126 ± 0.005) Mq. pc^^ for zq = 0. Considering spherically symmetric 



halo, its contribution to mass density is Qhaloi''') = (47rG) 



d-ir v^^iJ/dr, where 



''^halo is circular velocity of an object due to the gravity of the halo. One can use 
Eq. (7). We show another possibility, now. We have [vGaiaxyiR)]^ = l^diskiR)]'^ + 

{vhaloiR)? + [vbulgeiR)]'^, VGalaxy{R) = [A{R) - B{R)] R, V,Usk{R) = [A^skiR) - 



Bd^sk{R)] R- Now, we have Qhalo{Ro) = (47rG)-l il^ ^ {i?o ^«^Lo " 2(^^ - B'^)Ro + 

disk 

{A + B), [dvdtsk{R)/dR]R„ = - {Adisk + Bdisk) are used. We have g^aio = Qhaloi.Ro) 
= 3.94 X 10~^ Mq pc~^ . (As a comparison we can mention that the halo mass density 



2 (^disfe - Bli^k) Ro - 2«6«;ge [d-Vbulge {r)/dr]Rj,if the relations [dvGalaxy (R) I dR\ 
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according to the halo model of Daupholc et al. (1996) is 2.35-times greater than our 
result: ghalo{Ro) (Dauphole et al.) = 9.26 x 10^'^ Mq pc^^ ■) We can summarize 

[^Galaxy (R)] = [Vd^skiR)] + [VhaloiR)] + {Vbulge{R)] , 
VGalaxyiR) = [A{R) " B{R)]R , 

'VdiskiR) = [Adisk{R) - Bdisk{R)]R : 



Vbulgeir) = \jr 



dr 
GMb 



QhaloiR) = {'^'^G)-^[X{Galaxy) - X{disk) - XQmlge)] , 
X{Galaxy) = - [A{R) - B{R)] x [A{R) + iB{R)\ 

X{disk) = - [Aa^MR) - BdiskiR)] X [AdrskiR) + SBdiskiR)] 
Yfh I ^ - '"bulge (R) ( Vbulge {R) „ dVbulge (R) 

A {bulge) — I + z 



R \ R dR 
Y=(A + Bfa\ + {A-ZBfal + 

(^disfe + Bdiskj'^'^Aaisk + {Adisk - ^Bdisk)"^ f^B^isk ' 

Mb = 1.3955 X 10^° Mq , 
bb = 0.35 kpc , 
Adisk = Ad^sk{Ro) = (10.81 ± 1.70) km kpc'^ , 
Bdisk = BdiskiRo) = (- 9.97 ±0.34) km s'^ kpc'^ , 
Qdisk = (0.126 ± 0.005) Mq pc'^ , 
Qhalo = (0.004 ± 0.001) Mq pc-^ , 

Q = Qdisk + Qhalo > (22) 

if also the bulge model of Dauphole et al. (1996) is inserted. The contribution of the 
bulge to the total mass density is Qbuige I {Qdisk + Qhalo) = 1 x 10"'*. The main part 
of the bulge is concentrated inside a sphere of small radius and the value of Qhalo is 
not very sensitive to the form of the bulge potential, Vf^y^ig^ « ^jGM\Jr. 

The total value of the mass density, coming from the disk and the halo, is = 0.130 
Mq pc~^. The value obtained from the model of Dauphole et al. (1996), 0.143 Mq 
pc~^, is in 10% greater than our result. The value of g seems to be in a good agreement 
with the values (0.120 ± 0.008) Mq pc~^ and (0.138 ± 0.009) Mq pc~^ obtained from 
the results of Agekjan in 1962 and Agekjan and Ogorodnikov in 1974 (see Kulikovskij 
1985, pp. 158-159) under the assumption that average mass of the stars in the solar 
neighborhood equals to 1 Mq (it does not significantly differ from 1 Mq). 

As for {dQ/dR)iig, the model of Maoz (2007) yields [dg{R,z = 0)/dR]Rg disk 
= g'^.^i^ = = (-0.0360 ± 0.0037) Mq pc~^ kpc~^ for zq = 0. Our model yields 
ldQ{r)/dr]Rg halo = Qhalo = - 0.62 x 10"=^ Mq pc-^ kpc-\ 

Perhaps, the result represented by Eq. (21) can be a little improved using the fact 
that 



Q{Ro,z) = Qdisk (1 - M \z\) + Qhalo , 
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dQ{R, z) 



dR 



Qdisk 



(1 - M \z\) + Qhalo , 



Ro 
Qdisk 



(0.126 ± 0.005) Mq pc"-^ , 
Qhalo = (0.004 ± 0.001) Mq pc-^ , 
Qdisk = (- 0.0360 ± 0.0037) Mq pc~^ kpc 
Qhalo = - 0.0006 Mq pc-^ kpc-^ , 
u = 3.3 kpc~^ , 



(23) 



where the value of u follows from the model of the disk. 

Eqs. (20) and (23) are apphed to both bodies, the one with the coordinates {Xq, 
Yo, Zq) and the second with the coordinates {X, Y, Z). In application to the Solar 
System and the Oort cloud, we have the Sun and the comet. We have 



fZo 
d?Z 



- |4 7r 

- |4 TT 



G 



G 



Qd^ 



Qdisk 



k (i - 

(X - Xo) 



\ u \Zo\) 



2 « 1^0 

+ ^ 



+ Qhalo 



+ 2{A^ - 



B 



+ Qhalo 



+ 2U 



B 



Yo)} 



'■)} ^0, 

=)} Z 



[Qdisk ~ ^ ^ l^l) + S'halo] Z 



(24) 



Initial conditions for the case of the Sun are: Zq = 30 pc, dZo/dt = + 6 km s 
(information on the basic solar motion is used, see, e.g., Mihalas and McRae Routly 
1968, p. 101). Using X = Xq + ^, Y = Yq + r] and Z = Zq + C, Eq. (24) yields 



dt^ 



GMq 



c 



AnG^QC - edisk\u \\Zq + CI (^0 + C) - l^ol^o]} 
2{J^ - B^) C 



'^0 



4 TT G ( - ddisk \ u 1^0 + Cl) (^0 + i 



(25) 



if also acceleration from the body with the coordinates (Xq, Vbi -^o) (Sun) is taken 
into account. The last term of Eq. (25) assures that motion in a plane parallel to the 
plane of the galactic equator does not exist: the inclination of the comet relative to 
the galactic equatorial plane does not fulfill the relation i = 0, since Zq = docs not 
hold. The statement that "plane-parallel motions do not exist" is equivalent to the fact 
that initial conditions C = and dC,/dt = do not admit C = constant = during 
the motion. Thus, our approach differs from the approach of Bottlingor in 1924-1925 
reducing to Oort equations for the galactic equatorial plane (see Kulikovskij 1985, pp. 
90-91). 



6.2 Summary 

Summarizing the most relevant equations and results presented above, we obtain 



d^ 
df^ 



GMq 



£, + {A- B) [A + B + 2A cos (2 wot)] $ 
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- 2 A (A - B) sin (2 cuat) tj 

+ 2 {A- Bf (A - Ta Zl) Ro Zq cos (ujot) C , 

- 7? - 2 A (A - B) sin (2 wq*) ^ 

+ (A - B) [A + B - 2 A cos (2 ojQt)] V 

- 2 {A-Bf (A - Zg) i?o sin(a;ot) C , 




-^C^ + 2(^^-5^)] C 



— 4: G Q Zq [cos (wo^) ^ — sin (wot) r?] , 

- [4 TT G e + 2 - S^^] Zo , 



r 



A - S , 



(26) 



if the terms containing the quantity u are neglected and the relations Xq = Rq 
cos(— Wot) and Yq = Rq sin(— wot) are used. The numerical values of the relevant 
quantities are: 



A = (14.2 ± 0.5) km s"^ kpc~^ , 

B = {- 12.4 ±0.5) km s^^ kpc^^ , 
A = 0.124 fcpc"^ , 
r2 = 1.586 fcpc"* , 
g = (0.130 ± 0.005) Mq pc~^ , 

q' = (- 0.037 ± 0.004) Mq pc~^ kpc~^ . (27) 



Wc stress that g = g{Ro,z = 0) = 0.130 Mq pc ^, Rq = 8.0 kpc. The mass density 
in the solar neighborhood, 30 pc above the galactic equatorial plane, is g = g{Ro,z = 
0.03 kpc) = (0.117 ±0.005) Mq pcr^ , sec also Eq. (23). 

One must bear in mind that the values of A and B significantly differ from the 
lAU recommended values A = 15 km kpc~^ , B = — 10 km kpc~^. The 
following results should be for certain: 



if also the result for the mass density of Creze et al. (1998) and the lAU recommended 
values are considered. Unfortunately, there exist results which do not fit Eq. (28) (e.g., 
there are values of A and B which are out of the above presented intervals, see Clemens 



A £ (10.0, 15.0) km kpc~^ , 
B € {- 15.0, - 10.0) km s"^ fcpc"^ , 
g € (0.07,0.15) Mq pc~^ , 
g € (- 0.043, - 0.035) Mq pc~^ kpc 



-1 



(28) 



(1985)). 
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7 Galactic tide for Dauphole et al. (1996) model of GaleLxy 

Taking into account Dauphole et al. (1996) model of Galaxy, the relevant equations for 
galactic tide are: 



r 



C + (A-B) [A + B + 2Acos{2 toot)] ^ 



dt2 



+ {A -BY [r^D/^Jhl + Zl + r^D) Ro Zo cos(a;oi) C 
d^V __GMq ^ _ (^-B) sin(2 ujot) i 



r 

2 A {A- B) sin (2 coot) n 

\2 ' 



+ (A-B) [A + B - 2 A cos (2 wot)] r] 



- {A-Bf {FiD/^bl + Z^ + r2D) Ro Zo sin (wot) C 

$ = -^C- [4.G. + 2(^^-B^)]C 

— 4 n G g' Zo [cos (wot) ^ — sin (wot) r?] , 



[4-kG g + 2{A^ - B^)] Zo 



wo = - B , 
A = 14.25 km kpc~^ , 
B — — 13.89 km kpc~^ , 



Fid = 0.084 kpc^^ , 
= 0.008 fcpc"^ , 
g = 0.143 Mq pc~^ 



g = - 0.0425 Mq pc~^ kpc~^ , 

ba = 0.25 kpc . (29) 



8 Conclusion 

The paper treats the effect of the galactic tide on the relative motion of two close 
bodies moving in the galactic disk. Equations considering the effect of the galactic 
tide on motion of a body near the Sun are given by Eqs. (26). Values of the physical 
quantities are presented in Eqs. (23) and (27), or, the more rough values are given by 
Eqs. (28). Taking into account Dauphole et al. (1996) model of Galaxy, equations for 
galactic tide are also presented. 

The effect of the normal a— component is conventional. The x— and y— components 
of the acceleration comes not only from the x— and y— components of the position of 
a body, but also from its a— component of the position. Motion of the Sun normal to 
the galactic equator plays also an important role. 

The physical models presented hero will be used in a forthcoming paper for a more 
realistic study of orbital evolution of long-period comets. 
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